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Abstract 

This letter describes a direct search for pair produced magnetic monopoles in e+e" collisions. 
The analysis is based on 62.7 pb^^ of data collected with the OPAL detector at an average centre- 
of-mass energy of v^= 206.3 GeV. The monopole signal was assumed to be characterized by two 
back-to-back particles with an anomalously high ionization energy loss dE/dx in the tracking 
chambers. No evidence for production of monopoles was observed. Upper limits were obtained on 
the magnetic monopole pair-production cross-section (cr) in the mass range 45 GeV/c^ < uim < 
102 GeV/c^. The average limit is < 0.05 pb and is essentially independent of the magnetic 
monopole mass. The cross-section limit is derived at the 95% confidence level and is valid for 
spin-1/2 magnetic monopoles. 
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1 Introduction 

In 1931 Dirac linked the existence of magnetic monopoles (MMs) with the quantization of electric 
charge and postulated the relation between the elementary electric charge e of the electron and a 
basic magnetic charge g [1]: 

g=r^=ngD, n=l,2,... (1) 

where n is an unknown integer and g£) = hc/2e = 68. 5e is the unit Dirac magnetic charge (in the 
cgs system). If free quarks exist, Eq. [T]should be modified by replacing e with e/3, which effectively 
increases 5 by a factor of 3. There was no prediction for the monopole mass. A rough estimate, 
obtained assuming that the classical monopole radius is equal to the classical electron radius, yields 
fnM ~ ^^w-e/e^ ~ "n? ■ 4700me ~ • 2.4 GeV/c^. Since 1931, experimental searches for "classical 
Dirac" monopoles have been performed at nearly every new high-energy accelerator, employing a 
variety of direct and indirect methods [2]. By a classical (Dirac) monopole, we mean a particle 
without electric charge or hadronic interactions and with magnetic charge g satisfying the Dirac 
quantization condition (Eq. [1]) . 

Within the framework of Grand Unified Theories (GUT) of the strong and electroweak interac- 
tions, supermassive magnetic monopoles with masses m > 10^^ GeV/c^ could have been produced 
in the early Universe as intrinsically stable topological defects at a high energy phase transition that 
leaves an unbroken U(l) group [3]. At the present time, such monopoles could exist in the penetrat- 
ing cosmic radiation as "fossil" remnants of that transition. The detection of such particles would be 
one of the most spectacular confirmations of GUT predictions. The most stringent upper limits on 
an isotropic flux of GUT magnetic monopoles, assuming monopole masses tum > 10^^ GeV/c^, have 
been set by the MACRO experiment [4]. In some Grand Unified theories values of the monopole 
mass as low as 10^ GeV/c^ are allowed [5,6]. Although it is not yet possible to set direct limits at 
this mass scale, it is worthwhile to search in the accessible region at LEP energies. 

Searches for classical point-like monopoles have been performed mainly at high-energy accel- 
erators and in cosmic radiation experiments. Monopole searches have predominantly used either 
ionization or induction detection techniques. 
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Induction experiments measure the monopole magnetic charge and are independent of monopole 
mass and velocity. These experiments search for the induction of a persistent current within a 
superconducting loop [7]. Searches for magnetic monopoles using this method have been performed 
at the pp Tevatron collider assuming that produced MMs could stop, and be trapped and bound, 
in the matter surrounding the DO and CDF collision regions [8]. The same strategy has been used 
to search for magnetic monopoles produced in e"'"p collisions at HERA [9]. 

Ionization experiments rely on the large magnetic charge of monopoles to produce more ionization 
than an electrical charge travelling with the same velocity. For g = go and velocities f3 = (v/c) > 
10~^ a magnetic monopole behaves, in terms of ionization energy loss (dE/dx), like an equivalent 
electric charge with {ze)eq = QdP- The energy losses are thus very large 

{dE/dx)g = {gf3/e)\dE/dx)e (2) 

and Dirac magnetic monopoles would be easily distinguished from minimum ionizing electrically 
charged Standard Model (SM) particles [10-12]. Direct searches for magnetic monopoles using 
tracking devices were performed at pp and e"'"e~ colliders. Experiments at the Tevatron collider 
established cross section limits of about 2 x 10~ cm for MMs with mu < 850 GeV/c^ [13], while 
searches at LEP have excluded masses up to 45 GeV/c^ [14]. 

Indirect searches for classical monopoles have relied on the effects of virtual monopole/anti- 
monopole loops added to QED processes in pp and e'^e^ collisions [15,16]. Since the Standard 
Model Z^boson could couple to monopoles, assuming that the coupling between the TP and a MM 
pair is larger than for any lepton pair, the measurement of the TP decay width provides an indirect 
limit on MM production for niM < mz/^ [6,12]. 

This paper describes a direct search for MM pairs produced in e"^e~ — > MM(-f) reactions. 
The data were collected with the OPAL detector at the LEP accelerator at CERN. This search was 
primarily based on the dE/dx measurements in the tracking chambers. OPAL has a well established 
analysis to search for stable, long-lived, massive particles using the dE/dx signatures of individual 
charged particle tracks [17]. This analysis technique could not be used here because MMs are too 
heavily ionizing, resulting in charge saturation in the central jet chamber. Therefore, a new analysis 
method was developed based on hit information rather than reconstructed tracks. The analysis was 
sensitive to MMs with masses from 45 GeV/c^ up to the kinematic limit (about 103 GeV/c^). 

2 The OPAL Detector 

A description of the OPAL detector and its jet chamber can be found in reference [18]. Only a brief 
overview is given here. 

The OPAL detector operated at LEP between 1989 and 2000 and is now dismantled. The central 
detector comprised a system of tracking chambers, providing track reconstruction over 96% of the 
full solid angl^ inside a 0.435 T uniform magnetic field parallel to the beam axis. It consisted of 
a two-layer silicon microstrip vertex detector, a high-precision vertex drift chamber with axial and 
stereo wires, a large-volume jet chamber and a set of z-chambers measuring the track coordinates 
along the beam direction. 

^The OPAL right-handed coordinate system is defined such that the z-axis is in the direction of the electron beam, 
the X-axis points toward the centre of the LEP ring, and 6 and (j) are the polar and azimuthal angles, defined relative 
to the +z- and -l-x-axes, respectively. The radial coordinate is denoted by r. 
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The jet chamber (CJ) [19] is the most important detector for this analysis. The chamber, with 
a diameter of about 2m and a length of about 4m, was divided into 24 azimuthal sectors, each 
equipped with 159 sense wires. Up to 159 position and dE/dx measurements per track were thus 
possible. 

The CJ also provided the hardware trigger for monopole candidates. This trigger identified 
events with highly ionizing particles. Of the 159 sense wires of a sector, 36 wires were combined to 
define three groups with 12 wires each. One group was at an inner region, close to the e+e~collision 
axis. The other two groups were at central and outer regions. For each wire, hits from highly ionizing 
tracks were identified as those yielding an integrated signal above a threshold of 1250 counts in the 
Flash Analogues to Digital Converters (FADC). For comparison, a minimum ionizing particle yields 
about 200 FADC counts. Values slightly above 1000 FADC counts are typical for protons with a 
momentum of a few hundred MeV. If, within a group, more than 10 wires detected a high dE/dx 
hit, a decision bit was set. If this bit was set by all groups of a sector, the monopole trigger was 
fired. Using raw hit information of randomly triggered events, the monopole trigger was determined 
to have an efficiency greater than 99%. 

A lead-glass electromagnetic calorimeter located outside the magnet coil covered the full az- 
imuthal range with good hermeticity in the polar angle range of | cos 6\ < 0.984. The magnet return 
yoke was instrumented for hadron calorimctry covering the region | cos 6\ < 0.99 and was surrounded 
by four layers of muon chambers. Electromagnetic calorimeters close to the beam axis completed 
the geometrical acceptance down to 24 mrad on each side of the interaction point. These small-angle 
calorimeters were also used to measure the integrated luminosity by counting Bhabha events [20] . 

In order to trigger on the signal described in the introduction, only data collected when the 
monopole trigger was active were used. The data-set analysed here was recorded during the LEP2 
phase with an average centre-of-mass (cm.) energy of 206.3 GeV, and corresponded to a total 
integrated luminosity of 62.7 pb^^. 



3 Monte Carlo Simulation 

The signal reaction e''"e~ — > MM was simulated at ^/smc = 208 GeV for monopole masses (mju) 
of 45, 50, 55, 60, 65, 70, 75, 80, 85, 90, 95, 100, 101, 102, 103 and 104 GeV/c^ with Monte Carlo 
(MC) event samples. Each sample contained 1000 events. Small differences in the centre-of-mass 
energies between the OPAL data analysed {^JSmin= 203.6 GeV, ^ Smax= 207.0 GeV, for an average 
y/sdata= 206.3 GeV) and the signal MC samples (y^sjuc) have a negligible effect on the analysis. 
MM masses were scaled to the cm. energy with the equation: 



V'-'aata /o\ 
(3) 
SMC 

This scaling is valid since dE/dx (hence detection efficiency) is a linear function of mass. 

The very large value of the magnetic charge makes it impossible to use perturbative theory to 
calculate the MM production process. MMs were assumed to be spin 1/2 particles, produced from 
the e~''e~ initial state via annihilation into a virtual photon, which yields a monopolc-antimonopole 
pair with a uniform azimuthal distribution and with the typical fermion polar angle distribution 
oc (1 -hcos2 6'): 

e+e~ ^ 7* ^ MM (4) 
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Since magnetic charge cannot be simulated directly, MMs were simulated as heavy electrically 
charged fermions with an effective charge of {ze)eq = gnP (assuming n = 1). The specific ionization 
energy loss was computed according to Eq. [2j 

A magnetic monopole interacts with a magnetic field analogously to how an electron interacts 
with an electric field. The Lorentz force for a magnetic monopole carrying magnetic charge g is: 

F = g(B-vxE) (5) 



The GEANT3 [21] based OPAL detector simulation program [22] was used to simulate the 
behavior of the MMs in the OPAL detector. The routines to transport the particles through the 
magnetic field were modified such that over a given step the change in the momentum dp/dt of the 
monopole was obtained by solving analytically the differential equation: 



The solution describes the motion of a magnetic monopole in a uniform magnetic field. The 
trajectory is a parabola, accelerating in the direction of the magnetic field. In the plane perpendicular 
to the magnetic field the motion is along a straight line, in sharp contrast to electrically charged 
particles, which curve in this plane. 

We studied the effects of multiple scattering of the monopoles and the modelling of the electric 
field between the anode, cathode, and potential wires in CJ and found them to be negligible. 

A software emulation of the monopole trigger was used to study its efficiency. For the simulated 
monopole events, the trigger efficiency was found to be essentially 100%. 

The background was estimated using MC simulations of Standard Model processes, generated 
at -^5=206 GeV. Two-fermion events {Z^*/^* ffil) with / = e,fi,T,q) were simulated with 
KK2f [23]. For the two-photon background, the PYTHIA [24] and PHOJET [25] Monte Carlo 
generators were used for e'''e~gg final states and the Vermaseren [26] and BDK [27] generators for 
all e^e~/^Z~ final states. Four-fermion final states were simulated with grc4f [28], which takes into 
account interference between all diagrams. 

All generated signal and background events were processed through the full simulation of the 
OPAL detector. The same event analysis chain was applied to the simulated events and to the data. 



4 Data Analysis 



Magnetic monopoles would distinguish themselves by their anomalously high ionization energy loss 
in CJ and by the different plane of curvature of the trajectory in the magnetic field, compared to 
electrically charged particles. 

The large value of the specific energy loss (dE/dx) of a MM in the gas of the tracking detectors 
would induce a saturation in most of the wire hits. With the signals from both ends of the wire 
saturated, it is not possible to determine the z position from charge sharing. In this case the z 
position is set to zero by the reconstruction program. In the MC, most MM events are seen to exhibit 
a mean z-coordinate near zero, because of saturation effects. Rather than trying to reconstruct the 
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Description 


cut value 




_L(JLcLl Clicliji,c pel lllL yy^O J. 






Number of Tracks plus Clusters: 


< 18 


1 


The first hit wire: 


< 2 




Number of Tracks plus Clusters: 


< 4 


2 


Distance between the 2 sectors: 


> 8 


3 


Number of hits in overflow in HighSector: 


> 10 


4 


Z mean coordinate (CJ): 


< 50 cm 


5 


Charge per hit in the HighSector: 


> 3700 FADC counts 


6 


Charge per hit in the Seconds ector: 


> 3000 FADC counts 




Total charge per hit (CJ): 


> 2500 FADC counts 



Table 1: List of cuts applied to the data. 



MM tracks in 3 dimensions, events were examined for the characteristic MM pattern of ionisation 
in the sectors of the OPAL Jet Chamber. 

Pair-produced magnetic monopoles, e"''e~ MM(7), would be expected to be produced back 
to back with a characteristic pattern of hits in the jet chamber. This would have resulted in an 
azimuthal separation of about 12 sectors between the two sectors with the highest energy deposits, 
called HighSector and SecondSector, with little energy deposited elsewhere in the detector. 

Based on these considerations, events were rejected if the overall charge deposited on the sense 
wires normalised per hit was smaller than 1000 FADC counts, or if the total multiplicity of tracks 
plus clusters in the detector was greater than 18. The level of the FADC counts were based on gains 
and calibrations. We refer to these two cuts as the preselection, see Tabled! 

To reject some un-modelled events, further cuts were applied: the number of reconstructed tracks 
plus clusters had to be no more than 4 and the first wire hit in CJ had to be one of the first two wires 
(cut 1 in Table [T]) . Table [T] summarizes the other selection criteria. We required the HighSector and 
SecondSector to have an azimuthal separation of at least eight sectors (cut 2) and the number of 
hits in overflow in the HighSector to be larger than or equal to 10 (cut 3). Since the typical MM 
signature would exhibit a mean z-coordinate near zero, the average of the z coordinate in CJ was 
required to be less than 50 cm (cut 4). The deposited charge per hit in HighSector and SecondSector 
was required to be larger than 3700 FADC and 3000 FADC counts, respectively (cut 5 and cut 6) 
and the total charge per hit in all the CJ sectors to be larger than 2500 FADC counts (cut 6). 

The Standard Model background was dominated by Bhabha events and two-photon hadronic 
events, with a contribution from other two-photon events. The effect of the cuts on the samples at 
an average cm. energy of y^=206.3 GeV is shown in Table [2j After applying cut 1, there was poor 
agreement between data and MC (see Table [2]). This was because the data still contained remaining 
un-modelled backgrounds from beam-gas interactions, cosmic rays and detector noise. This un- 
modelled background was much reduced by the subsequent cuts since beam-gas interactions yield 
particles which mainly travel along the beam pipe and do not have the characteristic back-to-back 
pattern, and detector noise does not deposit large amounts of charge on the wires. The remaining 
difference of 15-20% between the number of events in the data and MC after cut 2 does not affect 
our results, as the signal is so separated from the the background that we can impose very hard cuts 
to remove all the background without affecting the detection efficiency. 

Fig. [1] shows the distribution of two of the main variables used by the analysis after cut 2: the 
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Numbe 


r of background events SM MC 






cut 


data 


Total SM MC 


bhabha 


2f 


qq 


27(e) 


27(m) 


27(t) 


w 


4f 


27(q) 


sig. cff.(%) 


1 


44491 


5707 


4231 


0.7 


0.6 


75.3 


2.2 


71.9 


1.9 


57 


1266 


91 


2 


2928 


2462 


1927 


0.1 


0.3 


6.0 


0.1 


27.5 


0.3 


14 


487 


91 


3 


2576 


2194 


1661 


0.1 


0.3 


5.4 


0.1 


27.4 


0.3 


12.8 


487 


91 


4 


1982 


1597 


1405 


0.0 


0.0 


0.4 


0.0 


6.9 


0.0 


6.4 


177 


91 


5 


2 


1.2 


0.5 


0.0 


0.0 


0.0 


0.0 


0.1 


0.0 


0.0 


0.6 


91 
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0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


91 



Table 2: The number of data and Monte Carlo events remaining after the cuts for anal- 
ysed data-set collected at y/s=206.3 GeV and for various MC SM background processes 
normalised to the integrated luminosity of the data (62.7 pb~^). The last column gives 
the efEciencies (in percent) for the magnetic monopole MC signal simulated in the mass 
region between 45 GeV/c? and 103 GeV/(?. 



charge per hit in the CJ sector HighSector and the average of the z-coordinate. The total number 
of data events at this stage is 2928 and the total number of the MC Standard Model events is 
2462 (Table [2]). Since the magnetic monopole behavior would be very different from any electrically 
charged SM particles, all the variables used by the analysis have a very well separated distribution 
for the MM signal and SM MC backgrounds. For this reason it can be seen from Table[2]that no MC 
background event survived the analysis cuts. Moreover the overall detection efficiency is very high 
(> 90%) for almost all MM masses. In Fig. [2] the detection efficiency for pair-produced magnetic 
monopoles at y/s = 206 GeV is shown as a function of mM- 



5 Estimates of Systematic Uncertainties 

The distributions of the variables in the data and SM MC have similar shapes. The differences in 
the mean values are quite small. The MC modelling of the dE/dx may introduce some systematic 
uncertainties. These were evaluated by displacing the cut value on a given variable x from the 
original position xq to a new position xq, to reproduce on the simulated events the effect of the cut 
on the real data, xq is defined by: 

xo = (xo - {x)^^J + {x)^^t^g (7) 

'J data 

where {x) 

bkg^ ^data and cjfefcg are the mean values and the standard deviations of the dis- 
tributions of the variable x for the data and the simulated background. These quantities were 
calculated from the x distributions of the events surviving the cuts on all the other variables used 
in the selection. It was verified that using the distribution of x at other stages of the selection leads 
to negligible changes in the values of this uncertainty. 

The procedure was repeated for the main variables used in the event selection (Table [T]): the 
number of overflows in HighSector, the Z mean coordinate in CJ and the charge per hit in HighSector 
and SecondSector. The difference between the reduced efficiency, due to the displacement of the cut, 
and that obtained with the nominal selection was taken as the systematic uncertainty due to the 
modelling of the variable under consideration. The relative systematic uncertainties in the signal 
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Quantity 


Systematic uncertainty (%) 


Number of overflows in HighSector 


0.0 - 0.2 


Z mean coordinate (CJ) 


0.2 - 0.4 


Charge per hit in HighSector 


0.3 - 4.7 


Charge per hit in SecondSector 


0.2 - 2.2 


Global systematic uncertainty 


0.4 - 5.2 


Signal MC statistics 


0.6 - 0.8 


Total 


0.7 - 5.3 



Table 3: Summary of systematic uncertainties for the signal efEciency of the various 
quantities used in the analysis. The range of results corresponds to the values obtained 
for the different MM masses. 



efficiency associated with the various quantities are reported in Table [3l The range comes from 
different values obtained for the different MM masses. 

At a given centre-of-mass energy the different systematic uncertainties were assumed to be 
independent, so that the total systematic uncertainty was calculated as the quadratic sum of the 
individual uncertainties. The global systematic uncertainty ranges between 0.4% and 5.2% (Table[3]). 

The MC statistical uncertainty, due to the limited number of signal events generated, has been 
computed using a binomial formula and is reported in Table [3l 

6 Results and Conclusions 

No magnetic monopole signal was found in this search. In Figure [3] the 95% CL upper limit on 
the production cross-section at an average cm. energy of ^/s = 206.3 GeV is shown as a function 
of the monopole mass. The average upper limit on the cross-section, computed using a frequentist 
approach, is 0.05 pb in the mass range 45 < mAf < 102 GeV/c^. This limit is essentially independent 
of the mass in this range. 

The computation of the cross-section is non-trivial. Nevertheless we expect the cross-section 
to be large. The cross-section for the pair production of Dirac Magnetic Monopoles computed 
assuming a naive tree-level coupling through an s-channel virtual photon, according to the effective 
charge {ze)eg = goP, is around 5 orders of magnitude larger than the upper limit obtained in this 
experiment [14]. In this model we can thus exclude classical MMs in the mass range 45-102 GeV/c^. 
This is a new excluded mass range for Dirac magnetic monopole searches in e"'"e~ interactions. 
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charge per hit in HighSector (FADC) average z-coordinate in CJ (cm) 



Figure 1: Global event properties after applying cut 2. The solid histogram is the generated Standard 
Model MC background, the filled histogram is the MC magnetic monopole signal and the points are 
the data, a) charge deposited in the HighSector in CJ (FADC); b) average z-coordinate of all hits 
in the jet chamber. Most MM events are seen to exhibit a mean z-coordinate near zero, because of 
the saturation effects mentioned in Sect. HI The arrows in each plot show the cuts applied in the 
analysis. 
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Figure 2: Monte Carlo estimate of the selection efficiency as a function of the monopole mass at 
^/s=206.3 GeV. The dotted line is the kinematic limit. 
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Figure 3: The model-independent 95% C.L. upper limit on the pair-production cross-section of 
magnetic monopoles in e^e~ collisions at LEP2 at y^=206.3 GeV (plotted vs MM mass). 
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